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Summary 

This paper discusses the Nastran Program Output from a structural analysts point of view, and 
describes the simple modifications which have been made to the program in order to improve it in 
certain areas. In particular, the convenience of the output for use in original design work is critically 
appraised and compared with the output from Astral.* It is shown that considerable hand calculation 
is necessary in order to extract useful load distribution data from the available Nastran output. For 
this reason, some effort has been directed towards providing additional force output for the Nastran 
shear panel element. 


Introduction 

The output from any finite element structural analysis program must be in a form which is 
immediately useful to the stress or design engineer. The most sophisticated program is not being used 
efficiently if lengthy hand calculations are necessary before useful loads data can be extracted from 
the computer output. It is for this reason that Astra! has been refined in a design environment over 
several years to provide a "Force” type output. This is the output that Grumman Engineers have 
found preferable to a stress output in most design situations where the structural model only 
approximates to the real structure. The output from Nastran is definitely not in a convenient form and 
this paper describes some of its shortcomings and the modifications made to the program in order to 
provide a force type output for the shear panel element. Similar modifications should also be made to 
the membrane and plate elements. 
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Shear panel width in. 

Shear panel depth in. 

Shear panel corner forces # 

Number of loading conditions 

Wing chord wise load distribution # /in. 

Shear panel diagonal forces # 

Shear flow # /in. 

Number of members in idealization 

Shear panel thickness in. 

Local coordinates in. 


| Angles defining geometry of shear panel (See Figure 73) 


NASTRAN OUTPUT: COMPARISON WITH ASTRAL 

In this section of the paper, Nastran Output is compared with Astral "M Print” for output format 
and form. 


* Automated Structural Analysis Program developed by Grumman 
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Output Format . Nastran output is grouped by loading conditions whereas Astral groups output by 
member or grid point number. Grouping by loading condition is far less convenient when many 
loading conditions may be present as is often the case in a fuselage analysis. A Iso, grouping by member 
number enables a maximum minimum search to be made by the program with a minimum of 
additional programming. This maximum minimum output for toads and deflections is an indispensable 
feature of Astral “M Print” output. 

As an example of the different grouping strategies, consider a hypothetical structure with r 
members and n loading conditions. Figure 1 compares the output formats for this structure. 

Output Form. Figure 2 indicates the main differences in output between Nastran and Astral for some 
of the most commonly used finite elements. In general, Nastran provides average stresses and scant 
toads data, whilst Astral provides a comprehensive force output on a ‘‘node force” basis. Nastran also 
provides a Margin of Safety output for some elements if an allowable stress is specified by the user. 
This margin is based upon average stresses in idealized members and assumes a uniaxial stress field. As 
such it is a ‘"misleading" output and should be disregarded in all but the most elementary truss type 
structure. Of more use, would be a stress ratio of the form: 


Stress Ratio 


Stress 

Allowable Stress 


Any stress output should be treated with caution, especially in final design work, since stresses 
output by Nastran are only correct if the idealized structure is identical to the real structure - rarely 
the case. In most cases, the Structural Analyst prefers to work with load distributions. 

Axial Force Output. Nastran outputs average forces in rod elements, diagonal forces in shear panels 
and no forces for the membrane element. Hence, the element node forces can not be calculated when 
using membranes, but can be calculated from rod and shear panel output at the expense of tedious 
hand calculations. The following procedure shows how node forces are obtained from rectangular 
shear panels and rods. Note that for the more general case of non rectangular panels, the required 
calculations are more tedious and an approximate graphical approach would be used in most cases. 


Assume a situation were a rod is bounded by two rectangular shear panels as shown in Figure 3: 


Nastran Output: 


Panel 1 Diagonal forces P] and P 2 

(#) 

Rod A verage Force P^ 

(#) 

Panel 2 Diagonal Forces P 4 and P$ 

(#) 
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Required: 


(#) 

(#/in) 


Node Forces at Rod Ends F-j and F 2 
Shear Flows Adjacent to Rod q-j and q 2 

Resolving diagonal forces along pane! edges, adding and dividing by length of rod: 


(P 1+ P 2 ) 



#/lN 


( p 4 + P 5) 



#/IN 


Summing corner forces and rod average force at each end of rod: 



These equations would be modified appropriately if more than two shear panels were adjacent to 
the rod. Typically, for a wing rib and beam intersection, four panels would be involved. 


H3-T Orbiter Wing Analysis 

This analysis has been chosen to illustrate how loads are extracted from both Nastran and Astral 
for use in design work. 

Briefly, the idealization consists of quadrilateral and triangular membranes for the upper and 
lower covers, and rods and shear panels for the ribs, beams (spars) and vertical studs (posts). A drawing 
of the idealization is shown in Figure 4. 

Main features of the idealization are as follows: 

130 Grid Points 
80 Quadrilateral Membranes 
14 Triangular Membranes 

107 Shear Panels 
71 Rods 

364 Degrees of Freedom 
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The Nastran Bulk Data Deck used for the analysis was generated automatically from the Astral 
Input Deck using an in house conversion program. Only the Control Decks were produced by hand. A 
comparison of computer running times for the two analyses is shown in Figure 5. 

Load distributions useful to the Structural Analyst and Designer are shown at sample wing 
sections, and the method of deriving them from the raw program output is detailed. The general 
agreement between the two analyses, as indicated by overall wing flexibility, is good as expected, since 
both use the stiffness method of finite element analysis and both have approximately equivalent 
elements. The difference in wing tip deflection for one typical loading condition amounted to less 
than .02%. 

Wing Cover Axial Loads 

Cover axial loads are most often presented at Grumman as a series of chordwise running load 
distributions(#/in.) along each rib tine. As an example, the distribution for the upper cover at rib 10 is 
determined in Figures 6 and 7 and shown plotted in Figure 8. The Astral node or cap loads are divided 
by the effective panel width. The Nastran average stresses in panels on the inboard and outboard sides 
of the rib must be converted to running load and then averaged. 

Note that both analyses require a finer grid in the area adjacent to the rear beam due to the 
discontinuity at grid points 91 and 92. Additional stringers between the beams which carry through to 
the airplane center tine are also needed, in order to predict the shear tag effect caused by the cover 
discontinuity across rib 1 1 . 

Rib and Beam Shear Flows 

AH shear flows are output directly by Astral “M Print”. Nastran output consists of maximum and 
average shear stresses which may be converted to shear flows by multiplying by the panel thickness, if, 
however, all four shear flows are required, then the procedure outlined previously for rectangular 
panels or the procedure below for extracting stud loads must be used. Figure 9 shows the extraction of 
maximum shear flows only, for rib 10. 

Stud Cap Forces 

Stud forces or cap forces are the most difficult to extract from Nastran output, in this example 
the cap forces for stud 85 at the intersection of rib 1 0 and beam 5 will be determined. None of the 
adjacent four shear panels are rectangular, therefore, the simplified equations of the Appendix 
may not be used, instead, an approximate graphical approach is taken. First the four shear panels 
are drawn to scale on Figure 10 with the diagonal forces and rod average forces as shown. Next, the 
components of the diagonal forces acting aiona the stud are scaled and added: 
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"167 "1 43 *1 39 

At the upper cover: Force = 3055 + 16292 x — 2514 x — : — +4393 x — 10911 x 

2.49 3.17 1.87 


1.78 

2.06 


= 6660# 


At the lower cover: Force = 3055 + 2626 x-^-- 15416 x -5-^+ 12199 x-^^-- 3960 x- — 6 - 

3.30 2.35 1.84 2.09 


= 1990# 

These stud forces are not in good agreement with the corresponding forces output directly by 
Astral which are 76/0# at the upper cover and 710# at the lower cover. This is probably due partly to 
the graphical approach used, and partly to the relative insensitivity of wing flexibility and deflection 
to stud loading. A Nastran wing solution with rigid studs, imposed by multipoint constraints, showed 
only a 4.6 % reduction in wing deflections. 

Program Modifications to Improve the Shear Panel Output 

This section of the paper describes in detail the changes made in Nastran Module SDR2 as a first 
step in providing an improved "force” type output. Perhaps the most inconvenient aspect of the 
present output is the use of shear panel diagonal forces, and it was decided to improve this output as 
quickly as possible. This has been accomplished by modifying the appropriate routines and 
subroutines in Module SDR2 and outputting the additional data directly from this module, without 
recourse to User Modules or changes in the Output File Processor. The additional output includes 
corner forces (#), shear flows(#/in), panel edge lengths and internal node numbers. This output in no 
way affects the regular Nastran printout which follows it: (See Appendix. ) 

In order to make these changes, the necessary equations were derived for resolving diagonal 
forces into corner forces and combining corner forces into shear flows. Next, a working knowledge of 
how Nastran interfaces "horizontally” between routines and " vertically ” between modules was 
obtained, mostly on a trial and error basis, since this information is not specifically given in the 
Programmers Manual. The layout chart for Module SDR2 shown in Figure 7 7 was then drawn and 
proved invaluable when modifying routines and deciding on the best data paths between routines. At 
first it was hoped to accomplish all data transfer by means of common blocks, however, it was soon 
realized that due to the looping that takes place within the module this was not possible. Instead, use 
was made of the data file ESTA when necessary. The necessary changes to module SDR2 are described 
in some detail below and summarized in Figure 12. (See Appendix.) 

The basic problem consists of bringing together the diagonal forces and local mean plane 
coordinates, performing the necessary calculations and printing the results. The diagonal forces are 
calculated in subroutine SPANL2 and the local coordinates in subroutine SPANL 7. A convenient place 
to calculate the forces and shear flows and print them is in routine SDR2E. The following list of 
modifications to the routines have been made. 



0 The local coordinates of the four corners of each panel are loaded into common block SDR2X5 
as SPOUT (8) in subroutine SPAN LI for transfer to SDR2B. These coordinates are X and Y 
values for corners 1 to 4 (Figure 13). 

0 In routine SDR2B, the local coordinates from common block SDR2X5 are transferred to Data 
File ESTA for the “vertical” transfer to SDR2E (Figure II). ESTA assembles a block of data for 
each shear panel element processed, whereas the common blocks are rewritten for each pass 
through a loop. 

0 The block data subroutine SDR2BD must be updated to redimension the enlarged ESTA file. The 
block size is increased from 25 to 33 words to include the coordinates. 

0 The shear panel diagonal forces PI 3 and P24 are loaded into common block SDR2X7 in SPANL2 
and transferred to SDR2E where they are available as BUFB(2) and BUFB(3) respectively. This 
does not require any additional coding. 

0 Data file ESTA is moved into common as ELESTA in routine SDR2E and element numbers and 
internal grid point numbers are already available as BUFB(I) and ELESTA(2-5) respectively. 
Hence, all the necessary data is assembled in SDR2E. 


Local Coordinates 
Node Number 
Diagonal Forces 
Member Number 
Subcase Number 


ELESTA (25-32) 
ELESTA (2-5) 
BUFB (2-3) 
BUFB (I) 
Z(ICC+1) 


0 In routine SDR2E the corner forces, shear flows and panel edge lengths are calculated. 
0 The calculated data is output with appropriate headings. 


Appendix A summarizes the equations used, shows a listing of the major SDR2E module coding 
changes and provides an example of the additional program output. 

Recommended Changes in Nastran Output 

The following changes in the form and format of Nastran Output would greatly ease the task of 
the Structural Analyst in design work. They would reduce the amount of hand calculation necessary, 
eliminate the likelihood of not designing to the most critical loading condition and increase the 
accuracy of most load distributions. 

Output Format 

0 The output should be grouped by member or grid point number, not by loading condition. 

0 The maximum and minimum loads for each member should be output if more than two loading 

conditions are present. 

0 Loads should be output in decimal form, not in exponential form, when of a suitable size. For 
example, a load of 5650# should be output as 5650.0 or 5650, not as 5.650000E 03 which is far 
more difficult to read correctly. 
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Output Form 

• The shear panel should have a shear flow (#/in.) output for all four sides, and the resolved corner 
forces should be used In a "Force ” type output including "kick ’’ loads due to warping, 

0 Rod and bar elements should have cap or grid point forces output, not average forces. These are 
determined using the above panel corner forces. 

0 The membrane element should have a "force" type output consisting of shear flows ( #/in .) along 
all four sides and cap forces along imaginary rod elements. 

0 The Margin of Safety output should be deleted. 

0 The length of all sides of all panels should be output, together with their direction cosines in the 
global coordinate system. 

0 A check on the free body equilibrium of the structure is essential output for any finite element 
analysis. This should be a summation of all forces and moments about the origin. 

Recommended Additions to Documentation 

0 Inclusion of Module Layouts in the Programmers Manual (see Figure II). 

0 Inclusion of full instructions for implementation of user modules in the programmers manual. 

Concluding Remarks 

The output from these two finite element analyses, Nastran and Comap Astral, differ in concept. 
Nastran in general outputs average stresses in elements and some forces which are not in an 
immediately usable form. Comap Astral provides no stress output (unless specifically requested) but 
gives a comprehensive force type output. Generally in all but the most elementary of structures, the 
force type output is preferred by structural analysts at Grumman. The Nastran Shear panel output is 
particularly difficult to use and the necessary modifications to the program to quickly provide a more 
suitable output have been described. At this time work is under way to provide a better output form 
and format for all the commonly used Nastran Elements. 

Reference 

AFFDL-TR-70-II8 An Automated Procedure for the Optimization of Pratical Aerospace Structures. 
Section 4.3 " Nodal Stress Method. " 
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Summary of Equations: 
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70 .4 80.400 .490. 490.490.490.500. 0000279 
.613 .610 .610.610. 610.610. 610. 610000260 

0000261 

0000262 

0000263 

0000264 

0000285 

0000266 

0000267 

ooooaee 

0000269. 


RES> 


P 24 = X8LFB <3 ) 

I !DO=8LF0 Cl> 

NODE 1 = (ELES7A (2)-l )/6 01 

NODE 2 ^CELESTA C3 )«» 1 ) /6 61 

NOOE3= CELESTA C4I-1 )/6Cl 
NODE 4 a CELESTA (5) -1 ) /6 Cl 

__ X2=XLESTA (2 7) 

X3»XLESTA (29) 

X4*XLE5TA C31) 

Y 2 = XLE £ T A <30 „ 

Y4=XLESTAC32) 

CA1 SY4/SQRT C >44X4 € Y4* Y4 ) 

_ SAlsX4/SQRTCX4*X46Y4»Y4) ' 

C A2 = Y 2/SC RT ( X 34 X 3 € Y3 ♦ Y3 ) 

£A2 = X2/£QRT C >3* X3 £Y34 Y3 ) 

CA3=Y2/SQRT CC X2-X3) 4426Y34Y3) 

SA3=C X2-X3) /SORTC CX2-X3)442CV34Y3) 
CA4«Y4/SQRT C C X2- >4) 4*2CY4 4 Y4 ) 

__ SA4s=C X2-X4) /SORT C CX2-X4 ) 4 426Y44Y4 > 

CA5 = C Y3-Y4) /SCRT C CX3-X4 ) 442&C Y3-Y4) 442) 

SA5 = ( >2- >4) /SORT ( ( X3- X4 ) 442&C Y3-Y4) 442) 

C C ALCULAT I CN CF CORNER FCRCES 

F 12=P134(SA24CA1-CA24£A1 )/C A1 
F 11 =P124CA2 /CAl 
F 22=P Z44CA4/CA3 

F 21=P244( SA44CA3-CA44SA3T/CA3 

F22=P134 C SA24CA3 6C A2 4 S A 3 ) / ( S A5 4C A3& C A54S A3) 
F3l=P124(<A54CA2-CA54SA2)/(SA54CA3& CAS4 S A3 ) 
F A2=P244< SA4 4C A5 6CA4 4SA5 )/ CSA54CA1-CA54S A1 ) 
F 4isP244C SA14CA4 €CA1 4SA4)/ CS/5 4CA1-CA54S A1 I 



APPENDIX - Concluded 


ISN 0239 
ISN 0240 
ISN 0241 
ISN 0242 
ISN 0243 
ISN 0244 
ISN 024* 
ISN 0246 
ISN 0247 
ISN 0246 
ISN 0 249 
ISN 0250 
ISN 0251 
ISN 0252 

ISN 0253 
ISN 0254* 
ISN 0255 
ISN 0 2 5 6 
ISN 0257 
ISN 0258 
ISN C259 
ISN 0260 
ISN 0261 
ISN 0262 
ISN 0263 
ISN 0 264 
ISN 0265 
ISN 0266 
ISN 0267 
ISN 0266 
ISN 0269 
ISN 0270 
ISN 0271 
ISN 0272 
ISN 0273 
ISN 0274 
ISN 0275 
ISN 0 276 
ISN 0277 
ISN Q27e 
ISN 0279 
ISN 0260 
ISN 0281 
ISN 0282 
ISN 0283 
ISN 0284 


ISN 0265 
ISN 0267 
ISN 0289 


C WRITE OUT DATA IN RCLT1NR SDR2E — FCLLOW5 LINE F4 1 = P24*< S A 1 * CA 4* • • • 
ALt = X2 

AL2 = SGRT( Y3*Y3£ (X2-X3>**2) 

AL3=SQRT( C Y3-T4 I * *2 6 t X3- X4 ) **2 ) 

AL4=SCPT { Y4 4Y4 6XA* X4 > 

S 1=<F 21-F 12 >/ALl 


S 2= (F 22— F 21 ) / A L2 
52 = (F4 1-F22 )/AL3 
£4 = CF42-F tl )/AL4 

WRITE C 6. 106) I IDO rNCDEl .NCDE2 *NCDE3 *N03E4, AL1# AL2# AL3# AL4 
10 6 FORMAT U X .1 6 .1 X.4 U10 «2>) .4Fl2*J| 

WRITE <6 • 1C7 ) F 1 1 «F2 1 «F 31 *F4 1 • S l * S2 • S3 »S4 
107 FOR MATC1CX*€F12*3) 

WR IT£(6«1C6) F12 % F22 tF32 ,F42 
ICE FORMAT! 10X.4FI2* 3 •//> 

C THIS COMPLETES DATA C ALCLL AT I C N AND OUTPUT 3/17/72 

GO TO 62 C 

450 K * 3 *" "" - 

GO TO 550 
4 6 C K* 0 

GO TO 510 ' •" 

47C Ks 3 

GO TO SIC 

4 60 K* l 

GO TO 5 30 

49C CALL SELAS2 

GO TD 620 
fOC K* 4 

. 510 CALL 5B5PL2IK) 

GO TO 620 

5 20 Ks 2 

■20 CALL «TQME2<K) 

GO TO 62 C 
•40 K * 4 

55C CALL STRCD2 (K) 

GO TO 62C ' 

560 CALL S5AR2 

GO TO 62 0 - 

57C AGAIN = .FALSE* 

CALL SCCNE2 ( SORC ) 

GC TO 62 C ; 

560 CAUL STRIR2 ITGRID) 

GO TO 620 

5 9C CALL STRAP2 <TGRIC ) 

GO TO 62 C 

6CC CALL STORC2 (TGRID ) 

GO TO 620 - __ _ 

610 GO TO 9CC 
C*** ** 

C CALL ELEMENT 2 TIMES FCR CCNFIEX VECTOR* IMAGINARY FIRST# REAL 
C CALL ELEMENT ROUTINE TWICE IF AXIC PROBLEM 
C ONCE FOR EACH OF THE 2 VECTORS IN CORE 

C***** 

620 IF( AXIC • A NO • MICVEC*NE*0 • AND* IP ART* EQ* 1 I GO TO 625 
IF! IPART * GE • KTYPE ) GC TC 615 

625 IVEC * MIOVEC 

C****4 

C F OR CONICAL SHELL CNLY 


0000290 
~ 0000291 

0000292 
0000293 
0000294 
0000295 
0000296 
0000297 
0000298 
0000299 
0000300 
0000301 
. 0000302 

0000303 
0000304 
00003C5 
0000306 
0000307 
0000308 
0000309 
0000310 
0000311 
0000312 
0000313 
0000314 
0000315 
0000316 
0000317 
oooo3ie 
0000319 
0000320 
0000321 
0000322 
SECOND0000323 
0000324 
0000325 
0000326 
0000327 
0000328 
QG00329 
0000330 
0000331 
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Output Format Comparison 


Nastran 

Loading Condition 

Member 


1 ! 

7 ■* r 


2 

1 - r 


n 

1 ■> r 

Astral 

Member 

Loading Condition 


1 

1 - n 


2 

1 - n 


r 

1 •* n 


FIG. 1 


Member Output Form Comparison 


Element Type 

Astral M Print Output 

Nastran Output 

Astral Bar 
Nastran Rod 

Node forces at each end of bar 

A verage force in bar 
A verage stress in bar 

Astral Beam 
Nastran Bar 

Node forces at each end of Beam. 
Moments, torque and shear at each 
end of beam 

A verage force in beam 
Moments, torque and shear 
as for Astral. AH respective 
stresses 

Shear Panel 

(Astral Type 6) 
(Nastran Shear) 

Shear flow (#/in.) along all four 
sides 

Length of each side 

Four diagonal forces on shear 
panel 

A verage and maximum shear 
stress #/in2 

Quadrilateral 
Membrane 
(Astral Type 8) 
(Nastran QDMEM) 

Node cap forces along each side 
Shear flow (#/in.) along each 
side 

Length of each side 



No force type output 
A verage stresses in local 
x and y directions and shear 
stress 

Principal stresses and maxi- 
mum shear stress 

Triangular 
Membrane 
(Astral Type 4) 
(Nastran TRMEM) 

Node cap forces only 

No force type output 
Stress output as for 
Quadrilateral membrane 


FIG. 2 
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SPAR 

9 


UPPER COVER ~ L.H. WING 

LOWER COVER NODES AND PANELS ARE 
NUMBERED ONE GREATER THAN SHOWN 

STUD NUMBERS ARE THE SAME AS THE 
UPPER GRID POINT NUMBERS 


SPAR 

2 



I 

I 

RIB 

1 


SCALE: APPROX, 1/100 


FIG. 4 ORBITER WING IDEALIZATION 


352 



Program Running Times - Comparison 


Time in Seconds 

Co map Astra 1 

Nastran 

CPU 

101.2 

214.3 

Job 

157.7 

241.3 


FIG. 5 
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Spanwise Load Intensity Distribution 


Astra! Results (M Print) 


Node 

Number 

Spanwise 
Load # 

Effective * 
Width in. 

N x 

#/in 

Distance 
From F. Bm in. 

79 

+1100 

35.63 

+31 

0 

81 

-4133 

70.66 

-59 

71.3 

83 

-69980 

91.57 

-764 

141.3 

85 

-130990 

93.66 

-1399 

254.4 

87 

-150400 

88.43 

-1701 

328.6 

89 

-118790 

86.62 

-1370 

431.3 

91 

-93280 

41.57 

-2240 

501.9 

93 

0 

6.26 

0 

514.4 


* Effective width is equal to the sum of half the panel widths on each side of 
the node. 


Fig. 6 


Nastran Results 


Panels 

tnbd of 

Panel 

N x 

Panels 

tnbd of 

Panel 

N X 

Average 

Distance 
Aft of 
Front Beam 

Rib 10 

F x * #/in 2 

t in 

#/in 

Rib 10 

F x *#/in 2 

t in 

#/in 

N x 

259 

-7680 

.025 

-192 

275 

-4210 

.0250 

-19 

-105 

35.6 

261 

-12990 

.025 

-325 

277 

-16510 

.0250 

-500 

-412 

106.3 

263 

-35350 

.025 

-884 

279 

-39820 

.0256 

-1134 

-1009 

197.8 

265 

-49500 

.0277 

-1371 

281 

-52060 

.0317 

-1732 

-1552 

291.5 

267 

-50000 

.0282 

-1410 

283 

-51800 

.0307 

-1646 

-1528 

379.9 

269 

-51400 

.0294 

-1511 

285 

-57820 

.0369 

-2371 

-1942 

466.6 


* F x is the average membrane stress in the panel, in the local coordinate system x direction. 


































Nastran and Astra! Maximum Shear Flows in Rib 10 


Panel 

Number 

Nastran Max . 
Shear Stress #/in. 

Panel 

Thickness in. 

Nastran Max. 
Shear Flow #/in. 

Astral Max. 
Shear Flow #/in. 

461 

11212 

.02 

224 

233 

463 

10473 

.02 

210 

244 

465 

2086 

.02 

42 

37 

467 

17297 

.02 

346 

327 

469 

9235 

.02 

185 

160 

471 

10915 

.02 

218 

289 

473 

19476 

.02 

389 

343 


FIG. 9 
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NASTRAN MODULE SDR2 


SDR2BD 

BLOCK 

DATA 


Fig. 11 


CASECC (IFP1) 
^XYCDB (IFP1) 


SDR2 


x 1» x 2> x 3 


SDR2AA 


, CASECC (IFP1) 


x 1> x 2» x 3 


SDR2A 


EDT (IFP) 
SIL(GPI) 
CSTM (GP1) 
H EST (TA1) 
GPTT (GP3) 


X 1 » x 2> x 3 

« ► 


SDR2B 


x 5> x 6 ^ 


tS 


SPAN LI 
ETC. 


(GP1) EQEXIN 
CASECC 
EIGR 

PG (SSG1)'\ 

y 


V 

'V. 


-svv CSTM (GP1)\\ 

v' 


x l» x 2> x 3 
«— 




SDR2C 


l— -31 , 

"TT — PUGV i 

1 1 PGPDT 

J|/,ESTA 

L ' 


UGV(SDRI) Vjj/ G PTTA 


CASECC EDT 
EIGR «. (IFP) ESTA \ 
PG GPTTA^' ' 


•x 


\ v x 2 , x 3 
4 - -» 


X4, x 7 , Xg 


SDR2D 


// 
ss 

a 

x 1> x 2» x 3 

< *> 


% ESTA 
^ 


X4, x 7 , Xg 


\N 

oesi n V s 
oefi N \^ 

V 

w!// 


SDR2E 


x 4> x 7 , x 8 

4 — ► 


SPANL2 

ETC. 


/> \\ N corner 

,/ OES1 N s FORCES 
*✓*' OEFI ** 

V 
M* * 

x- 


ADDITIONAL DATA FOR 
CORNER FORCE MODIFICATIONS 
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Summary of Changes to SDR2 Module to Provide Shear Panel Corner Force Output 


Routine 

Location 

Change or Inclusion 

SDR2BD 

Line 88 
"Data Elem 1 ” 

Increase "ESTA ” block for "shear” element by number of words 
being added (8 in this case) column 40 ... 33 ( was 25) 

SDR2E 

Line 36 

Add: "Dimension XLESTA (100), XBUFB (200)” 

These are ‘real’ array equivalents of integer arrays 'ELESTA ’ 
and BUFB 

SDR2E 

Line 96 

Add: " Equivalence (ELESTA (1), XLESTA (1) ), (SBUFB (1), 
BUFB (1))” 

SDR2E 

Line 111 

Add Headings for each subcase. This output will appear at the 
beginning of each subcase. 

SDR2E 

Line 289 

Add: a) calculations for shear flows, etc. 
b) output of shear flows, etc. 

SPAN LI 

Line 69 
" common ” 

In common " (SDR2X5 f” and SPOUT (8) (to hold corner 
coordinates) 

Continuation card 3: "3, SPOUT (8), yyyyy (94)” 

SPAN LI 

Line 191 

Add: SPOUT (1)„„ SPOUT (8) -0.0, 0.0, X2, 0.0, X3, Y3, X4, 
Y4 


Fig . 12 
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